Background cosmological dynamics in f{R) gravity and observational constraints 
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In this paper, we carry out a study of viable cosmological models in /(7?)-gravity at the background 
level. We use observable parameters like and 7 to form autonomous system of equations and show 
that the models under consideration exhibit two different regimes in their time evolution, namely, 
a phantom phase followed by a quintessence like behavior. We employ statefinder parameters to 
emphasize a characteristic discriminative signature of these models. 
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I. INTRODUCTION 

One of the most challenging problems of modern cos- 
mology today is associated with the attempts to under- 
standing the late time acceleration of Universe which is 
supported by cosmological observations of complimen- 
tary nature such as Supernovae la [l| , Cosmic Microwave 
Background anisotropics 0, Large Scale Structure for- 
mation jSj, baryon oscillations [j] and weak lensing [5J 
observations. Many theoretical approaches have been 
employed to explain the phenomenon of late time cosmic 
acceleration. The standard lore assumes the presence of 
an exotic fluid known as dark energy. The simplest dark 
energy model based upon cosmological constant dubbed 
ACDM model suffers from extreme fine-tuning and coin- 
cidence problems [6|. Scalar fields minimally coupled to 
gravity, called quintessence, with generic features might 
allow to alleviate these problems Many other possi- 
bilities have been proposed, including a scalar field with 
a non-standard kinetic term (/c-essence) |,8hl6l|, or sim- 
ply an arbitrary barotropic fluid with a pre-determined 
form for p{p), such as the Chaplygin gas and its various 
generalizations (T7l - [23 |. 

As an alternative to dark energy, the large scale mod- 
ifications of gravity could account for the current ac- 
celeration of universe. We know that gravity is modi- 
fied at short distance and there is no guarantee that it 
would not suffer any correction at large scales where it is 
never verified directly. Large scale modifications might 
arise from extra dimensional effects or can be inspired 
by fundamental theories. They can also be motivated by 
phenomenological considerations such as f{R) theories of 
gravity (see [25] for a recent review). However, any large 
scale modification of gravity should reconcile with local 
physics constraints and should have potential of being 
distinguished from cosmological constant. 

Most of the f(R) gravity models proposed in the lit- 
erature either ruled out by cosmological constraints im- 
posed by the history [2?'] or fail to meet the local gravity 
constraints |42i]. The viable f(R) models can be distin- 
guished from the ACDM by studying the the evolution of 
the growth of matter density perturbations [1^, [2^, H^l • 

In this paper, we explore the possibility of discriminat- 
ing the ACDM from viable models of f{R) at the back- 
ground level (see also Ref.[29| on the related theme). 



II. f{R) COSMOLOGY 

In what follows, it would be convenient to us to write 
f(R) gravity action in the form 
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167rG 



d4xV^[i? + e(i?)]+5,„(g^„*™) (1) 



where G is the bare gravitational constant, e{R) is a 
function of the curvature scalar R only and Sm is a func- 
tional of some matter fields and metric g^jy. 
In case of fiat homogenous and isotropic universe 



(2) 



the action ([T|) gives rise to the following evolution equa- 
tions (we use the unit, SttG ~ 1) 
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The equations (I3l4p can be rewritten with the defini- 
tion of the density of dark energy 
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-2H = Pr, 
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where wde = Pde/ Pde- 

In order to study the late time evolution in the f{R) 
models under consideration, we introduce the set of vari- 
able (51, 7, i?) (SI is the ratio of the dark energy density 
and the critical density) and 7 = \ + wde 



n' = 3Sl(l-f7)(l-7) 
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where a prime indicates differentiation with respect to 
In a. One can see this an autonomous system of equa- 
tions involving the observable cosmological parameters 
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Figure 1: The left panel (a) shows the time evolution of model (1) for different values of the parameters (n, A) in the phase 
space {Q,-y). For all parameters, the model is close to ACDM in the past followed by a phantom phase and the amplitude of 
this phase depends on the parameters of the model. For large values of (n, A), the model in indistinguishable from the ACDM 
model. The transition occurs between the phantom and non-pantom phases around fl < 0.5 (for the range of parameters 
studied) before the matter-DE equivalence. The epoch of this transition depends on the parameters of the model. For large 
values of these parameters, the transition shifts towards the large values of Q corresponding to small values of redshift. 
The right panel (b) shows the time evolution of the model (2) for various values of A. For all models, the de Sitter point 
(n,7) — (1,0) is an attractor. 



like 7 and Q. 



III. COSMOLOGICAL EVOLUTION 



In what follows, we shall be interested in the f(R) 
models which contain a standard matter phase [27] to 
be compatible with the early universe physics (BBN) 
(e(i?) — > Constant in the past) and give rise the de- 
Sitter attractor at late times. Consistency also demands 
that these models should satisfy the stability criteria 
((1 + e,H > 0); {e"(i?) > ofisl) and be consistent with 
the local gravity constraints (4^] . 

Bearing in mind the aforesaid, we consider the follow- 
ing models: 

• (A) e{x) = -Ai?c + 1) [H, 

• (B) eix) = -Ai?e (1 - (1 + 

• (C) e(x) = ~XRc (1 - e-^) [H, 

• (D) e{x) ^ -~XRc tanh(a;) [U' 

where x — R/Rc and Rc of the order of the observed 
cosmological constant. 

Investigation reveals that models (A), (B) and (C), 
(D) are cosmologically duplicate. Thus in the analysis to 
follow, we shall focus on models (A), (C) and would refer 
to them as model (1) and model (2) respectively. 



We have investigated models (1) & (2) numerically. To 
do this, we have assumed that initially the models are 
closed to ACDM with 7i = and is negligibly small. 
We find that for all viable models belonging to these two 
classes, we have the same evolution in the (7, il) plane. 
Fig. (HI shows that for both models, we have a phantom 
phase in the past (small D, = Que) giving rise to violation 
of weak energy condition. As shown in Fig.([T]), both the 
models suffer a transition across the phantom divide in 
the past (57 < 0.7). For the models (1), with smaller 
values of n, the models deviate from ACDM whereas for 
models (2), the same is true for parameter A. 



IV. STATEFINDER ANALYSIS 

As demonstrated in (37| . it is possible to discriminate 
different models of dark energy from each other using the 
statefinder parameters (r, s) [36| , 
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Figure 2: The left panel (a) shows the time evolution of the statefinder pair {s,r} for the model (2). We have the same 
evolution for various values of the free parameter A. All models outline loops around the ACDM model. The upper left part of 
this plane corresponds to the past evolution of the system where the model shows a phantom phase and the lower right part is 
the non-pantom phase. 

The right panel (b) shows the time evolution of the pair {q, r}. The solid line is the time evolution of the ACDM model which 
divides the surface into 2 planes. The upper part is the phantom evolution of the model while the lower half is the non-pantom 
phase of the model. For all models the de-Sitter (dS) point, (s,r) = (0, 1) or equivalently {q,r) — ( — 1, 1), is an attractor. 



where r is the jerk parameter and s is a function of the 
jerk and the decelerating parameter (g) 



(12) 



3(g-l/2) 



The statefinder parameters are a natural next step be- 
yond the Hubble function H = a/ a. By adding more 
derivative of the scale factor, Sahni et al. constructed 
geometrical parameters (defined using the metric only) 
which can discriminate various models. This method can 
be used to distinguish our models from ACDM model, 
characterized by {r,s) — (1,0). 

For both types of models (1) and (2), we found the 
same evolution in the (s,r) and {q,r) planes. The mod- 
els under consideration are close to the ACDM model in 
the past (see Figl2]) which appears clearly in the (r, s) 
plane and shows that the system is close to the crit- 
ical point (r = l,s = 0). The violation of the weak 
energy condition, which corresponds to the early evolu- 
tion of the system, defined a loop in the upper left of 
the (r, s)-plane. In fact for a simple model like a power 
law evolution of the scale factor a{t) ~ t^^^'* , we have 
r — {1 — 37)(1 — 37/2) and 5 = 7. Then a violation of 
the weak energy condition (7 < 0) is characterized by 
r > 1 and s < 0. The system crosses the phantom line 



(FigH]) given by 7 = which corresponds to the ACDM 
point in the (r, s)-plane. An another loop is exhibited by 
the model for the quintessence evolution of the system 
which corresponds to r < 1 and s > in the simplest 
power law model. 

It must be emphasized that the aforesaid features of cos- 
mological evolution appear for all the viable models in 
f{R) that we have studied. These models generically 
exhibit two different dynamical regimes, a phantom evo- 
lution in the past followed by quintessence like phase at 
late time. 



V. OBSERVATIONAL CONSTRAINTS 

We constrain the free parameters of the models stud- 
ied by using Supernovae data and the BAO data. We 
used the compiled Constitution set [sT'l of 397 type la 
supernovae for which the is defined by 
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Figure 3: 68%, 95% and 99.9% confidence intervals for the model (1). In the left panel (a) we imposed Q,m,o ~ 0.2 and 
f2m,o = 0.3 for the right panel (b). 
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Figure 4: 68%, 95% and 99.9% confidence intervals for the 
model (2). The local gravity constraints and the stability of 
the de-Sitter phase impose A > 1. 
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^J'th^ = Slog {dL{zi)) + ^iQ + -rlog 



Gcs{zt) 



where /xq = 25 + 5 log 



■H- 



Geff(z = 0) 

is marginalised [H, H 



(14) 



and dL is the luminosity-distance. 

The addition of the last term in (|14l) takes into account 
a varying gravitational constant [35]. We will not 
include this term in the numerical analysis. In fact it is 
negligible in viable /(i?)-gravities models because of the 
thin-shell effect. 

We also used the BAO distance ratio Dy{z = 
0.35)/L»i,(z 0.2) = 1.736 ± 0.065 [13, where 



D,{z) 







H(z) 


Wo H{z')) _ 



1/3 



(15) 



In case of model (1) we fixed the value of fi^.o (Fig. 
[3]) and find that the model is strongly sensitive to this 
parameter. For r^m.o = 0.2, the small values of the 
parameters (n, A) are preferred. This is the range of 
the scalar regime (26j which is crucially different from 
ACDM model. While the model is totally unconstrained 
for ilm^o = 0.3. 

We observe that the model 2 is unconstrained by the 
data (FigH]) . The density of matter today is constrained 
around the concordance value and A appears like a free 
parameter. 

We also use information criteria (IC) to assess the 
strength of models. These statistics favors models that 
give a good fit with fewer parameters. We use the 
Bayesian information criterion (BIG) and Akaike infor- 
mation criterion (AIC) to select the best fit models. The 
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Figure 5: AIC(left panel)and BIC(right panel) versus^lmfi for the standard case of cosmological constant and for the model 2 



AIC and BIG are defined as 

AIC = -2lnL + 2k (16) 

BIG = -2lnL + klnN (17) 

Where L is the maximum likelihood, k is the number 
of parameters and N is the number of data used in the 
fit. For gaussian errors, — —2lnL, we plot the best 
fit values of the AIC and BIG as the function of flm.o 
for the standard model based on cosmological constant 
and model 2 respectively (Fig. 5) (For model (1) this has 
been already done in [44|. We can see from the figure 
that while the cosmological constant gives slightly better 
fit, when larger or smaller values of flm.o Siie considered 
the AIG and BIG tests shows model 2 is slightly favoured 
over AGDM. . For model 2, AAIC = 2 and ABIC = 6. 

VI. CONCLUSION 



that in the early time the models are close to AGDM. 
As a generic feature, these models exhibit two distinct 
regimes in the time-evolution of the system. In the past, 
the models violate the weak energy condition a la the 
phantom phase followed by a quintessence like behavior 
at the present epoch. A simple analysis at small red- 
shift can lead to wrong conclusion that the quintessence 
model fits the data perfectly. Fitting our models with 
Snia data as well as with the BAO, we see that except 
for small Qm where there is strong bound on n for models 
(1), there is no significant constraints on both the models 
from cosmological observations. 

We emphasize that a comprehensive analysis should be 
conducted independently in two different ranges of the 
redshift, i.e., after and before the equivalence between 
dark-energy and matter. In our opinion, this could be a 
smoking gun for models which with different evolution- 
ionary phases of dark energy. 



In this paper, we have examined the cosmological dy- 
namics of two different classes of viable /(i?)-gravity 
models. To begin with, we have formed an autonomous 
system of equations involving observable cosmological pa- 
rameters like 7 and f2 together with the Ricci scalar R. 
This is interesting set of equations and can give interest- 
ing results upon phase-plane analysis. 

We then use this system equations to study the cos- 
mological behaviour of the models. For this we assume 
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